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The equal-time pairing correlation function of the two-dimensional t-J model on a square 
lattice is studied using a high-temperature expansion method. The sum of the pairing correla- 
tion, its spatial dependence, and the correlation length are obtained as functions of temperature 
down to T ~ 0.2t. By comparison of single-particle contributions in the correlation functions, 
we find an effective attractive interaction between quasi-particles in d ;c 2_ H 2-wave pairings. It 
is shown that d-wave correlation grows rapidly at low temperatures for 0.5 < n < 0.9, with n 
being the electron density. The temperature for this growth is roughly scaled by J/2. This is in 
sharp contrast to the Hubbard model in a weak or intermediate coupling region, where there 
is no numerical evidence of superconductivity. 
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Since the discovery of high- T c superconductors, 1 great 
efforts have been devoted to explaining the mechanism 
of their superconductivity. It is now believed that two 
dimensionality and strong electron correlation are im- 
portant for high-T c superconductivity. Although many 
studies have been carried out using the Hubbard model 
and the t-J model on a two-dimensional square lattice, 
numerical studies using these models do not always sup- 
port superconductivity. In the Hubbard model, quantum 
Monte Carlo simulations show that there is no enhance- 
ment of pairing correlation, 2-5 although RPA, 6 fluctua- 
tion exchange (FLEX) approximations 7 or renormaliza- 
tion group approaches 8 claim d 2 ,2_ ?/ 2-wave superconduc- 
tivity. It is considered that the lattice size or electron 
correlation (U/t < 4) in the numerical calculation is too 
small to deal with superconductivity which has a very low 
energy scale. 9 Although several efforts to detect super- 
conductivity have been made, 10 it is far from conclusive 
whether superconductivity occurs in the two-dimensional 
Hubbard model numerically. 

In the t-J model, on the other hand, mean-field the- 
ory indicates d x 2_ y 2-w&ve superconductivity and super- 
conductivity is supported by several numerical calcula- 
tions. 11-19 However, these studies are restricted to small 
system sizes or a variational theory in the ground state. 
In exact diagonalization studies of small clusters, 12 ' 13 it 
is very difficult to discuss the long-range order. In vari- 
ational calculations, 14-18 on the other hand, it is known 
that the variational energy is mainly determined from 
short-range correlations, and thus it is difficult to con- 
firm the existence of superconducting long-range order. 
Therefore, some exact numerical calculations are highly 
necessary. Since quantum Monte Carlo simulations have 
not been successful in the t-J model, the only avail- 
able exact calculation is a high-temperature expansion 
method, which has no negative sign problem and can be 
used to treat large systems exactly. Using this method, 
we can study the properties at finite temperatures which 
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give complementary information for ground-state stud- 
ies. 

Using high-temperature expansions, thermodynamic 
quantities, and two-point correlation functions, such as 
free energy and spin-spin correlation functions, have been 
studied. 20-22 However, superconducting pairing correla- 
tions have not been obtained because calculation of these 
quantities, which are essentially four-point correlations, 
requires a greater number of graphs. Recently, Pryadko et 
al. have obtained a high-temperature series of pairing 
susceptibilities through the ninth order in inverse tem- 
perature and claimed that the pairing susceptibility does 
not grow at low temperatures for J/t < 1. However, the 
order of the expansion is too low and the error bars are 
too large to make a definite conclusion. 

In this study, we investigate the high-temperature ex- 
pansion of equal-time pairing correlation functions up 
to the twelfth order in inverse temperatures. In fact, se- 
ries analysis of a pairing correlation function is a dif- 
ficult problem since pairing correlation functions decay 
very rapidly for long distances as discussed below. To 
overcome this difficulty, we develop a method of observ- 
ing superconductivity on the basis of correlation lengths. 
In this analysis, we find that it is possible to study su- 
perconductivity more accurately. Using Pade approxi- 
mations, we successfully extrapolate the series down to 
T ~ 0.2£ and find that the pairing correlation for d x i_ y i- 
wave symmetry grows rapidly at low temperatures for 
J/t ~ 0.4 and 0.5 < n < 0.9. By this calculation, wc 
can explicitly show the development of superconducting 
correlations in the t-J model at finite temperatures. 

The t-J model is defined as 



P 



j 2^ s» • Sj , 



i) 



where summations are over the nearest-neighbor pairs on 
a square lattice and the projection operator P eliminates 
doubly occupied sites. For this model, we calculate the 
high-temperature scries of the equal-time pairing corre- 
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lation function given as 

P a (i,j) = (Al(r J )A Q (r,) + A Q (r,)At(r,)). (2) 

Here, a denotes the symmetry of the pairing function, 
and the order parameter A Q (rj) is defined by 



A Q (ri) = y^/q(Q(ci T Ci-m ±Ci+ifCii), 



(3) 



with / a (Z) being the form factor of the pairing correlation 
and Yi being the coordinate at site i. Considering the 
nearest neighbor pairings, the form factors in the square 
lattice are given by 



/s(0 = 7^,o(<^,i +Si v + ^„,o(<^,i +S tx _i) 



/d(0 = g^.o^.i+^v- 1 )- ^,0(^,1 +^x,-i) 



(4) 



(5) 



/ P (Z) = -=[^,0(^,1-^,-1)], 



(6) 



pairing corre- 
-wave and 



for a singlet case, and 
1 

vt L 

for a triplet case. Each denotes the 
lation for an extended s-wave, c? a ,2_ y 2-wave and p- 
wave, respectively. The average in eq. (2) is taken 
over the grand canonical ensemble; that is, (O) = 
Tr0e -/3(w-, 1 ]v e ) / / Tre -/3(H- / ,Jv e )_ We ca i cu i ate p a (ij) U p 

to the twelfth order in inverse temperature, (3, using the 
finite cluster method. 24 Since we consider Cooper pairs 
of nearest-neighbor sites, graphs appearing in the cal- 
culation have up to fourteen edges. In each graph, a 
pair of edges represents Cooper pairs at site i and j, 
and the rest of the edges represent the corresponding 
nearest-neighbor bonds in the Hamiltonian. After cal- 
culating P a (i,j) as a function of 0, A = e^, J and t, 
we eliminate A by solving the high-temperature series of 
electron density, n, given by 

dfl 2A 

<9yU 



+ 



1 + 2A 

and obtain P a (i,j) as a function of n. 25 The obtained 
series is too long to show here, but it is available from 
the authors upon request. 

For extrapolating a high temperature series to low 
temperatures, we use a Pade approximation. To improve 
the convergence, we also calculate the series using an al- 
ternative choice of the variable w(/3) = tanh(//3), where 
/ is a parameter tuned to optimize the convergence. 21 
For small values of /, the series does not differ much 
from that for (3. Note that when [3 — > 00, w becomes fi- 
nite, which avoids the unphysical divergences or zeros of 
a Pade approximation at (3 = 00. However, one cannot 
extrapolate to values of j3 much larger than 1// since 
tanh(//3) almost saturates. With this in mind, we gener- 
ate various Pade approximations using 0.01 < / < 0.2. 
The extrapolated values and the error bars of all the fol- 
lowing figures are calculated as averages and deviations 
of various Pade approximations, respectively. 



s-wave 




Fig. 1. Sum of equal-time pairing correlation functions, S a (T), 
for a = s,p, and d at n = 0.8 and J/t = 0.4. 



Let us first discuss the sum of the equal-time pairing 
correlation function (q = component), S a (T), defined 
as 



N 



J)- 



(7) 



Figure 1 shows S a (T) (a = s,p, and d-wave) at J/t = 0.4 
for n = 0.8. If superconductivity is realized at T = 0, 
S a (T) diverges as T — > together with pairing suscep- 
tibility. It can be seen in Fig. 1 that p-wave pairing cor- 
relation does not show an enhancement at low temper- 
atures. This is the expected behavior since the effect of 
J, which favors singlet pairing, is large near half filling. 
In contrast, s- and d-wave pairing correlations grow at 
low temperatures (T < 0.5i). In particular, s-wave pair- 
ing increases rapidly, which at first glance suggests an 
s-wave superconductivity. However, as we show shortly, 
we find that this s-wave enhancement is due to a local 
singlet formation. 

The pairing correlations as functions of distance, 
P a (r), given by 



P a (r) 



|r 4 -p^|=r 



(8) 



-Tj\=r 

with a — s or d are shown in Fig. 2. Here, the distance 
is defined as |r| = \r x \ + \r y \. In agreement with the 
behavior of Sd, Pd(f) grows with decreasing temperature. 
Figure 2(a) shows that Pd(r) with r > 2 will reach values 
of 0.01 or larger as T — > 0, indicating that the long- 
range pairing correlation exists at T = 0. In contrast, the 
s-wave pairing correlation becomes very small (or even 
negative in some parameters) for long distances (r > 3). 
This means that the enhancement of S s in Fig. 1 comes 
only from the short-range pairing correlation and that S s 
does not diverge as T — > 0. These results are consistent 
with the ground state studies. 

Although the result of P a (r) shows an enhancement 
of the long-range d-wave pairing correlation, the extrap- 
olation of P a (r) can be done only at relatively high tem- 
peratures (T/t > 1) for large r values, where the high- 
temperature series is rather short. For example, the se- 
ries of P a {r) at r — 5 starts from the order of /3 8 . Thus, 
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Fig. 2. Distance dependence of pairing correlation function, 
P a (r), at J/t = 0.4, n = 0.75 for (i) a = d and (ii) a = s. In the 
case of d-wave, Pd(r) with r = 1 ~ 5 grow with decreasing tem- 
peratures though absolute values are small for long distances. In 
contrast, P s (r) with r > 3 are not found in this figure since P s (r) 
at r = 4 or 5 is negative in this case. 



an alternative method is desired to clarify the growth 
of long-range correlation at low temperatures more ac- 
curately. For this purpose, we calculate the correlation 
length of pairing correlations as follows. Let us consider 
the Fourier transform of P a (i,j), 
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Fig. 3. Correlation lengths for s- and ci-wave pairing at J/t = 
0.4 and n = 0.8. For comparison, the correlation length of d- 
wave pairing for the free-electron system, £?, and the correlation 
length of single-particle contribution (see the text), are also 
shown. 



and does not grow at low temperatures, which also agrees 
with the behavior of P s (r) . 

Although the enhancement of £d is promising for su- 
perconducting long-range order at T = 0, we must be 
careful with our conclusion since the pairing correla- 
tion contains single-particle contributions. For example, 
the correlation length of cZ-wave pairing symmetry in 
the free-electron system, f° r the same electron den- 
sity (n = 0.8) is shown in Fig. 3. Although £|] is rel- 
atively large and grows as T — > 0, there should be no 
instability of superconductivity. This enhancement of £° 
comes from the single-particle part: for free systems, 
(c^Cjicl^) = (citCjT)( c J'i c lj.) ho-ldS; and the single- 
particle part (cijcj|) provides the enhancement. 

To extract an intrinsic pairing correlation, we use a 
method employed in the Hubbard model. 26 The single- 
particle contribution in P a (i,j) can be obtained from 
high-temperature series of (c^cj^) and (cjj,cj^). (They 
correspond to disconnected diagrams.) By applying the 

same extrapolation scheme used for £ d , we calculate £ d 

,t 



^a(q) = jj^2Pa(i,j)exp(ici- (r, - r.,)). (9) f rom^ (c^) (c^cj^). Tncn > comparing £ d with £ d , we 



When the correlation of A a develops, P a (q) shows a peak 
at q = 0. Then, the correlation length of A Q is defined 
from the width of its peak as 

d 2 

r^(q) 



£ = - 



1 



2rfP Q (q) 9q 2 
1 



9 =0 



E 



r i | 2 <At(0)A a (r < ) + A a (0)At(r i )), 

(10) 

with d being a dimension. We find that the high- 
temperature series of this quantity converges much bet- 
ter than P a (r). Figure 3 shows £ s and £ d at n = 0.8. 
It is apparent that £ d shows a large enhancement for 
T < 0.1-0.2t, which is consistent with the behavior of Sd 
and Pd{r). In contrast, £ s almost saturates at T ~ 0.4i 



can discuss the effective interactions among pairs. The 
obtained results using the t-J model are shown in Fig. 
3. It is observed that £ d is strongly suppressed compared 
with which indicates that electrons have difficulty 
moving with strong correlations. Furthermore, £ d , which 
is very similar to £ d at high temperatures, deviates from 
£ d and starts to grow at low temperatures. This means 
that there is an effective attractive interaction in the 
d-wave pairings. The growth of correlation length itself 
means that the pairing correlation becomes long-ranged 
with lowering of the temperature. In fact, the behavior 
of £ d is apparently different from that of the s-wave case, 
£ s , where pairing correlations are short-ranged. Although 
the divergence of £ d as T — > is not perfectly proved, the 
behavior of £ d at T < 0.3t strongly suggests supercon- 
ductivity in the ground state. 

The doping rate dependences of £ d and £ d are shown 
in Fig. 4 for some values of T/t. Apparently, £ d shows 
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Fig. 4. £d and ^ as functions of n at T/t ■ 
(from top to bottom). 



: 0.4,0.8,1.2, and 1.6 



enhancement over a wide region. However, in the low- 
density region, £<j and £<j are very similar since the ef- 
fects of J and the Gutzwillcr projection, which makes a 
difference between and £d, are small. We can see that 
£d grows faster than £d in the range 0.5 < n < 0.9. Thus, 
superconductivity is expected in this region. This doping 
dependence is consistent with that observed in previous 
ground state studies. 11, 14, 17 

Here, we describe the J dependence of pairing corre- 
lations. It is found that the behavior of £<j for J = <~ 1 
does not change qualitatively. Even at J = 0, grows 
with lowering temperatures, though it does not show 
rapid growth at T > 0.1 <~ 0.2t. With increasing J, the 
growth of £d starts at higher temperatures. It is roughly 
scaled by J/2 at n = 0.8. 

To summarize, we have studied the equal-time pairing 
correlation function of the two-dimensional t-J model 
on a square lattice. It is found that the pairing corre- 
lation grows with decreasing temperature. A prominent 
enhancement of correlation length is detected in a wide 
region, which occurs at T < 0.2t. This growth of the 
pairing correlation length indicates the superconducting 



ground state in the t-J model. This is in sharp contrast 
to the Hubbard model in a weak or intermediate cou- 
pling region. Although the pairing susceptibility in the 
Hubbard model shows an effective attractive interaction 
in the (i-wave channel, 26 pairing correlations are short- 
ranged and do not show scaling with lattice size. 3 ' 4 



10 

ii 

12 
13 

14 
15 
16 
17 
18 

19 
20 

21 

22 

23 

24 

25 
26 



J. Bednorz and K. Muller: Z. Phys. B 64 (1986) 189. 
M. Imada and Y. Hatsugai: J. Phys. Soc. Jpn. 58 (1989) 3752. 
N. Furukawa and M. Imada: J. Phys. Soc. Jpn. 61 (1992) 3331. 
A. Morco: Phys. Rev. B 45 (1992) 5059. 

S. Zhang, J. Carlson and J. E. Gubcrnatis: Phys. Rev. Lett. 
78 (1997) 4486. 

D. J. Scalapino: Phys. Rep. 250 (1995) 329. 

N. E. Bickers, D. J. Scalapino and S. R. White: Phys. Rev. 
Lett. 62 (1989) 961. 

C. J. Halboth and W. Metzner: Phys. Rev. B 61 (2000) 7364. 
H. Yokoyama, Y. Tanaka, M. Ogata and H. Tsuchiura: J. 
Phys. Soc. Jpn. 73 (2004) 1119. 

K. Kuroki and H. Aoki: Phys. Rev. B 56 (1997) R14287. 
M. Kohno: Phys. Rev. B 55 (1997) 1435. 

E. Dagotto and J. Ricra: Phys. Rev. Lett. 70 (1993) 682. 

E. Dagotto, J. Riera, Y. C. Chen, A. Moreo, A. Nazarenko, F. 

Alcaraz and F. Ortolani: Phys. Rev. B 49 (1994) 3548. 

H. Yokoyama and H. Shiba: J. Phys. Soc. Jpn. 57 (1988) 2482. 

C. Gros: Ann. Phys. (N.Y.) 189 (1989) 53. 

H. Yokoyama and M. Ogata: J. Phys. Soc. Jpn. 65 (1996) 3615. 

A. Himcda and M. Ogata: Phys. Rev. B 60 (1999) R9935. 

S. Sorella, G. B. Martins, F. Becca, C. Gazza, L. Capriotti, 

A. Parola and E. Dagotto: Phys. Rev. Lett. 88 (2002) 117002. 

E. Dagotto: Rev. Mod. Phys. 66 (1994) 763. 

W. O. Putikka, M. U. Luchini and T. M. Rice: Phys. Rev. Lett. 

68 (1992) 538. 

R. R. P. Singh and R. L. Glenister: Phys. Rev. B 46 (1992) 
11871. 

W. O. Putikka, R. L. Glenister, R. R. P. Singh and H. Tsunet- 
sugu: Phys. Rev. Lett. 73 (1994) 170. 

L. P. Pryadko, S. A. Kivelson and O. Zachar: Phys. Rev. Lett. 
92 (2004) 067002. 

C. Domb and M. S. Green: Phase Transitions and Critical 
Phenomena, vol. 3 (Academic, London, 1974). 
K. Kubo and M. Tada: Prog. Theor. Phys. 71 (1984) 479. 
S. R. White, D. J. Scalapino, R. L. Sugar, N. E. Bickers and 
R. T. Scalettar: Phys. Rev. B 39 (1989) R839. 



